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Summary. — Entanglement is not only the most intriguing feature of quantum 
mechanics, but also a key resource in quantum information science. Entanglement is 
central to many quantum communication protocols, including dense coding, telepor- 
tation and quantum protocols for cryptography. For quantum algorithms, multipar- 
tite (many-qubit) entanglement is necessary to achieve an exponential speedup over 
classical computation. The entanglement content of random pure quantum states is 
almost maximal; such states find applications in various quantum information pro- 
tocols. The preparation of a random state or, equivalently, the implementation of a 
random unitary operator, requires a number of elementary one- and two-qubit gates 
that is exponential in the number Uq of qubits, thus becoming rapidly unfeasible 
when increasing Uq. On the other hand, pseudo-random states approximating to 
the desired accuracy the entanglement properties of true random states may be gen- 
erated efficiently, that is, polynomially in Uq. In particular, quantum chaotic maps 
are efficient generators of multipartite entanglement among the qubits, close to that 
expected for random states. This review discusses several aspects of the relationship 
between entanglement, randomness and chaos. In particular, I will focus on the fol- 
lowing items: (i) the robustness of the entanglement generated by quantum chaotic 
maps when taking into account the unavoidable noise sources affecting a quantum 
computer; (ii) the detection of the entanglement of high-dimensional (mixtures of) 
random states, an issue also related to the question of the emergence of classical- 
ity in coarse grained quantum chaotic dynamics; (iii) the decoherence induced by 
the coupling of a system to a chaotic environment, that is, by the entanglement 
established between the system and the environment. 

PACS 03.67.Mn - Entanglement production, characterization, and manipulation. 

PACS 03 . 67 . Lx - Quantum computation. 

PACS 03 . 67 . -a - Quantum information. 

PACS 05 . 45 . Mt - Quantum chaos; semiclassical methods. 

PACS 03.65.Yz - Decoherence; open systems; quantum statistical methods. 
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2 G. BENENTI 

1. Introduction 

Entanglement [1, 2, 3], arguably the most spectacular and counterintuitive manifes- 
tation of quantum mechanics, is observed in composite quantum systems. It signifies the 
existence of non-local correlations between measurements performed on particles that 
have interacted in the past, but now are located arbitrarily far away. We say that a two- 
particle state IV') is entangled, or non-separable, if it cannot be written as a simple tensor 
product |A:i)|A:2) = |fci)(g)|fc2) of two states which describe the first and the second subsys- 
tems, respectively, but only as a superposition of such states: {ip) = X]fc k Cfcifea |fci)|fe)- 
When two systems are entangled, it is not possible to assign them individual state vectors. 

The intriguing non-classical properties of entangled states were clearly illustrated by 
Einstein, Podolsky and Rosen (EPR) in 1935 [4]. These authors showed that quantum 
theory leads to a contradiction, provided that we accept (i) the reality principle: If we 
can predict with certainty the value of a physical quantity, then this value has physical 
reality, independently of our observation; (^) (ii) the locality principle: If two systems are 
causally disconnected, the result of any measurement performed on one system cannot 
influence the result of a measurement performed on the second system. (^) The EPR 
conclusion was that quantum mechanics is an incomplete theory. The suggestion was that 
measurement is in reality a deterministic process, which merely appears probabilistic 
since some degrees of freedom (hidden variables) are not precisely known. Of course, 
according to the standard interpretation of quantum mechanics there is no contradiction, 
since the wave function is not seen as a physical object, but just as a mathematical tool, 
useful to predict probabilities for the outcome of experiments. 

The debate on the physical reality of quantum systems became the subject of ex- 
perimental investigation after the formulation, in 1964, of Bell's inequalities [5]. These 
inequalities are obtained assuming the principles of realism and locality. Since it is 
possible to devise situations in which quantum mechanics predicts a violation of these 
inequalities, any experimental observation of such a violation excludes the possibility of 
a local and realistic description of natural phenomena. In short. Bell showed that the 
principles of realism and locality lead to experimentally testable inequality relations in 
disagreement with the predictions of quantum mechanics. 

Many experiments have been performed in order to check Bell's inequalities; the most 
famous involved EPR pairs of photons and was performed by Aspect and coworkers in 
1982 [6]. This experiment displayed an unambiguous violation of a Bell's inequality 
by tens of standard deviations and an excellent agreement with quantum mechanics. 
More recently, other experiments have come closer to the requirements of the ideal EPR 
scheme and again impressive agreement with the predictions of quantum mechanics has 
always been found. Nonetheless, there is no general consensus as to whether or not these 
experiments may be considered conclusive, owing to the limited efficiency of detectors. 
If, for the sake of argument, we assume that the present results will not be contradicted 
by future experiments with high-efficiency detectors, we must conclude that Nature does 



C') For example, if a system's wave function lip) is an eigenstate of an operator A, namely, 
A\xf}) — a\^), then the value a of the observable A is, using the EPR language, an element of 
physical reality. 

(^) Following the theory of relativity, we say that two measurement events are causally discon- 
nected if (Aa:)^ > c^(Af)^, where Ax and At are the space and time separations of the two 
events in some inertial reference frame and c is the speed of light (the two events take place at 
space-time coordinates {xi,ti) and {x2,t2), respectively, and Ax — X2 — xi, At = t2 ~ ti). 
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not experimentally support the EPR point of view. In summary, the World is not locally 
realistic. 

I should stress that there is more to learn from Bell's inequalities and Aspect's ex- 
periments than merely a consistency test of quantum mechanics. These profound results 
show us that entanglement is a fundamentally new resource, beyond the realm of clas- 
sical physics, and that it is possible to experimentally manipulate entangled states. A 
major goal of quantum information science [7, 8] is to exploit this resource to perform 
computation and communication tasks beyond classical capabilities. 

Entanglement is central to many quantum communication protocols, including quan- 
tum dense coding [9], which permits transmission of two bits of classical information 
through the manipulation of only one of two entangled qubits, and quantum teleporta- 
tion [10], which allows the transfer of the state of one quantum system to another over 
an arbitrary distance. Moreover, entanglement is a tool for secure communication [11]. 
Finally, in the field of quantum computation entanglement allows algorithms exponen- 
tially faster than any known classical computation [12]. For any quantum algorithm 
operating on pure states, the presence of multipartite (many-qubit) entanglement is nec- 
essary to achieve an exponential speedup over classical computation [13]. Therefore the 
ability to control high-dimensional entangled states is one of the basic requirements for 
constructing quantum computers. 

Random numbers are important in classical computation, as probabilistic algorithms 
can be far more efficient than deterministic ones in solving many problems [14]. Random- 
ness may also be useful in quantum computation. Random pure states of dimension N 
are drawn from the uniform (Haar) measure on pure states (^) The entanglement content 
of random pure quantum states is almost maximal [16, 17, 18] and such states find appli- 
cations in various quantum protocols, like superdense coding of quantum states [19, 18], 
remote state preparation [20], and the construction of efficient data- hiding schemes [21]. 
Moreover, it has been argued that random evolutions may be used to characterize the 
main aspects of noise sources affecting a quantum processor [22] . Finally, random states 
may form the basis for a statistical theory of entanglement. While it is very difficult to 
characterize the entanglement properties of a many-qubit state, a simplified theory of 
entanglement might be possible for random states [18]. 

The preparation of a random state or, equivalently, the implementation of a random 
unitary operator mapping a fiducial n^-qubit initial state, say ]0) = ]0 . . . 0) = ]0) (8) • ■ • (8) 
]0), onto a typical (random) state, requires a number of elementary one- and two-qubit 
gates exponential in the number Uq of qubits, thus becoming rapidly unfeasible when 
increasing Uq. On the other hand, pseudo-random states approximating to the desired 



(^) The Haar measure on the unitary group U{N) is the unique measure on pure A'^- level states 
invariant under unitary transformations. It is a uniform, unbiased measure on pure states. For 
a single qubit (A^ = 2), it can be simply visualized as a uniform distribution on the Block 
sphere [7, 8]. The generic state of a qubit may be written as 

(1) |i/)) = cos-|0) +e''^sin-|l), < 61 < tt, < < 27r, 

where the states of the computational basis {|0), jl)} are eigenstates of the Pauli operator CTz. The 
qubit's state can be represented by a point on a sphere of unit radius, called the Bloch sphere. 
This sphere is parametrized by the angles 9, (j) and can be embedded in a three-dimensional 
space of Cartesian coordinates (a:: — cos <j) sin 0,y — sin (j> sin 0, z = cos 9) . 
We also point out that ensembles of random mixed states are reviewed in [15]. 
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accuracy the entanglement properties of true random states may be generated efficiently, 
that is, polynomially in Uq [22, 23, 24, 25, 26, 27]. In a sense, pseudo-random states 
play in quantum information protocols a role analogous to pseudo-random numbers in 
classical information theory. 

Random states can be efficiently approximated by means of random one- and two- 
qubit unitaries [22, 23, 24, 25, 26, 27, 28] or by deterministic dynamical systems (maps) 
in the regime of quantum chaos [29, 30, 31, 32, 33, 34, 35, 24]. These maps are known to 
exhibit certain statistical properties of random matrices [29, 30] and are efficient gener- 
ators of multipartite entanglement among the qubits, close to that expected for random 
states [35, 24]. Note that in this case deterministic instead of random one- and two-qubit 
gates are implemented, the required randomness being provided by deterministic chaotic 
dynamics. A related crucial question, which I shall discuss in this review, is whether the 
generated entanglement is robust when taking into account unavoidable noise sources 
affecting a quantum computer. That is, decoherencc or imperfections in the quantum 
computer hardware [7, 36], that in general turn pure states into mixtures, with a corre- 
sponding loss of quantum coherence and entanglement content. 

This paper reviews previous work concerning several aspects of the relationship be- 
tween entanglement, randomness and chaos. In particular, I will focus on the following 
items: (i) the robustness of the entanglement generated by quantum chaotic maps when 
taking into account the unavoidable noise sources affecting a quantum computer (Sec. 7); 
(ii) the detection of the entanglement of high-dimensional (mixtures of) random states, 
an issue also related to the question of the emergence of classicality in coarse grained 
quantum chaotic dynamics (Sec. 8); (iii) the decoherencc induced by the coupling of a 
system to a chaotic environment, that is, by the entanglement established between the 
system and the environment (Sec. 9). In order to make this paper accessible also to 
readers without a background in quantum information science and/or in quantum chaos, 
basic concepts and tools concerning bipartite and multipartite entanglement, random and 
pseudo-random quantum states and quantum chaos maps are discussed in the remaining 
sections and appendixes. 

2. Bipartite entanglement 

2T. The von Neumann entropy. - In this section, we show that for pure states \ip) a 
good measure of bipartite entanglement exists: the von Neumann entropy of the reduced 
density matrices. Given a state described by the density matrix p, its von Neumann 
entropy is defined as 

(2) S{p)^-TT{p\ogp). 

Note that, here as in the rest of the paper, all logarithms are base 2 unless otherwise 
indicated. 

First of all, a few definitions are needed: 

Entanglement cost: Let us assume that two communicating parties, Alice and Bob, 
share many Einstein-Podolsky-Rosen (EPR) pairs (*), say 

(4) [</-+) = i= (100) + 111)), 



(^) A basis of (maximally) entangled states for the two-qubit Hilbert space is provided by the 
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and that they wish to prepare a large number n of copies of a given bipartite pure state 
1^), using only local operations and classical communication. If we call fcmin the minimum 
number of EPR pairs necessary to accomplish this task, we define the entanglement cost 
as the limiting ratio k^in/n, for n —* oo. 

Distillable entanglement: Let us consider the reverse process; that is, Alice and Bob 
share a large number n of copies of a pure state |^) and they wish to concentrate entan- 
glement, again using only local operations supplemented by classical communication. If 
^max denotes the maximum number of EPR pairs that can be obtained in this manner, 
we define the distillable entanglement as the ratio k'^^^/n in the limit n — ^ oo. 

It is clear that k'^^^ < fcmin- Otherwise, we could employ local operations and classical 
communication to create entanglement, which is a non-local, purely quantum resource 
(it would be sufficient to prepare n states |'0) from fcmin EPR pairs and then distill 
^max > ^min EPR states). Furthermore, it is possible to show that, asymptotically in 
n, the entanglement cost and the distillable entanglement coincide and that the ratios 
kmin/n and k'^^^/n are given by the reduced single-qubit von Neumann entropies. Indeed, 
we have 

(5) lim ^ = lim ^ = S{pa) = S{pb), 

n^oo n n^oo fl 

where S{pa) and S{pb) are the von Neumann entropies of the reduced density matrices 
PA = TrB(|i/))(^|) and pB = Tr^(|i/))(^|), respectively. Therefore, the process that 
changes n copies of \'ip) into fc copies of 10"*") is asymptotically reversible. Moreover, it 
is possible to show that it is faithful; namely, the change takes place with unit fidelity 
when n —> oo.(^) The proof of this result can be found in [38]. We can therefore quantify 
the entanglement of a bipartite pure state |'0) as 

(9) EabM) ^ SipA) = SipB). 

It ranges from for a separable state to 1 for maximally entangled two-qubit states (the 



four Bell states (EPR pairs) 



(3) |0±)^_L(^|oo)±|ii)), |v;±) = -L(^|oi)±|io)). 

(^) The fidelity F provides a measure of the distance between two, generally mixed, quantum 
states p and a: 

(6) F(p, a) = (^Tr\/pi/2crpi/2^ ^ . 
The fidelity of a pure state {ip) and an arbitrary state a is given by 

(7) F(|V),a) = {vi|a|V), 

which is the square root of the overlap between |i/>) and a. Finally, the fidelity of two pure 
quantum states |?/>i) and |'!/;2) is defined by 

(8) F(|^l),!^2)) = l(V'l!^2)|'. 

We have < F < 1, with F = 1 when \ipi} coincides with |-!/)2) and F = when |V>i) and |'!/)2) 
are orthogonal. For further discussions on this quantity see, e.g., [7, 8]. The average fidelity 
between random states is studied in [371. 
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EPR states). Hence, it is common practice to say that the entanglement of an EPR pair 
is 1 ebit. 

2'2. The Schmidt decomposition. ~ The fact that S{pa) ~ S{pb) is easily derived 
from the Schmidt decomposition. 

The Schmidt decomposition theorem: Given a pure state \'ip) € Ti. = Ha 'S) Ti-B of a 
bipartite quantum system, there exist orthonormal states {|?)a} for T-La and {|i')B} for 
T-Lb such that 

k 

(10) IV) = ^ ^|i)A|0i3 - ^^|1)A|1')B + • • • + VP~k\k)A\k')B, 

i=\ 

with Pi positive real numbers satisfying the condition X]i=i Pi = 1 (for a proof of this 
theorem see, e.g., [7, 8]). 

It is important to stress that the states {|i)yi} and {|«')b} depend on the particular 
state \'^) that we wish to expand. The reduced density matrices pA = ^iViV) A a{t\ 
and pB = X^iPiNOs _b(«'| have the same non-zero eigenvalues. Their number is also 
the number k of terms in the Schmidt decomposition (10) and is known as the Schmidt 
number (or the Schmidt rank) of the state {ip). A separable pure state, which by definition 
can be written as 

(11) m = \^)a\Ob, 

has Schmidt number equal to one. Thus, we have the following entanglement criterion: 
a bipartite pure state is entangled if and only if its Schmidt number is greater than one. 
For instance, the Schmidt number of the EPR state (4) is 2. 
It is clear from the Schmidt decomposition (10) that 

(12) S{pa) = S{pb) = -Y.P^\ogp,■ 

i 

If N , Na and Nb denote the dimensions of the Hilbert spaces 7i, TLa and TLb, with 
N = NaNb and Na< Nb, we have 

(13) < EABim = S{pa) = S{pb) < log Na. 

A maximally entangled state of two subsystems has Na equally weighted terms in its 
Schmidt decomposition and therefore its entanglement content is logTV^ ebits. For in- 
stance, the EPR state (4) is a maximally entangled two-qubit state. Note that a maxi- 
mally entangled state \4') leads to a maximally mixed state pA- 

2'3. The purity. - The purity of state described by the density matrix p is defined as 

(14) P{p) = Tr(p2). 
We have 

(15) P{pa)^P{pb)^Y.pI 
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The purity is much easier to investigate analytically than the von Neumann entropy. 
Moreover, it provides the first non-trivial term in a Taylor series expansion of the von 
Neumann entropy about its maximum. (®) The purity ranges from 1/Na for maximally 
entangled states to 1 for separable states. One can also consider the participation ratio 
^ = v=^— 3", which is the inverse of the purity. This quantity is bounded between 1 and 
Na and is close to 1 if a single term dominates the Schmidt decomposition (10), whereas 
^ = Na if all terms in the decomposition have the same weight (pi — ■ ■ ■ — p^^ — 
1/Na)- The participation ratio ^ represents the effective number of terms in the Schmidt 
decomposition. 

A natural extension of the discussion of this section is to consider bipartite mixed 
states, p = J2iPi\''P){''P\ iJ2iPi — 1)' instead of pure states. However, mixed-state entan- 
glement is not as well understood as pure-state bipartite entanglement and is the focus 
of ongoing research (for a review, see, e.g., Refs. [1, 2, 3]). 

By definition, a (generally mixed) state is said to be separable if it can be prepared by 
two parties (Alice and Bob) in a "classical" manner; that is, by means of local operations 
and classical communication. This means that Alice and Bob agree over the phone on the 
local preparation of the two subsystems A and B. Therefore, a mixed state is separable 
if and only if it can be written as 

(17) PAB ^'^PkPAk® PBk, with pfe > and ^pfe = 1, 



where pAk and pBk are density matrices for the two subsystems. A separable system 
always satisfies Bell's inequalities; that is, it only contains classical correlations. Given a 
density matrix pab , it is in general a non-trivial task to prove whether a decomposition as 
in (17) exists or not [2, 3]. We therefore need separability criteria that arc easier to test. 
Two useful tools for the detection of entanglement, the Peres criterion and entanglement 
witnesses, are reviewed in Appendix A. 

3. Entanglement of random states 

A simple argument helps understanding why the bipartite entanglement content of a 
pure random state ji/i) is almost maximal. In a given basis {\i)} the density matrix for 
the state \tfj) — X]i<^il*) i^ written as follows: 

(18) p,, = (i|^)(^|j) = c,c*, 

where c^ — {i\ip) are the components of the state \tp) in the {\i)} basis. In the case of 
a random state the components are uniformly distributed, with amplitudes Ci w 1/yN 
and random phases. Here N is the Hilbert space dimension and the value 1/\/N of 
the amplitudes ensures that the wave vector ji/i) is normalized. The density matrix can 



(^) If we write pi = -^^, with e^ ^ 1 and ^^ et = 0, we obtain 
(16) S{pa) ~ logiVA - [NA/{2\n2)]P(pA). 
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therefore be written as 

(19) fl« diag (—,—,...,— I +0, 

where D, is a N x N zero diagonal matrix with random complex matrix elements of 
amplitude « 1/iV. Suppose now that we partition the Hilbert space of the system into 
two parts, A and B, with dimensions Na and Nb, where NaNb — N. Without loss of 
generality, we take the first subsystem. A, to be the one with the not larger dimension: 
Na < Nb- The reduced density matrix pA is defined as follows: 



(20) pa = TtbP = Y^ c,^iB 4^ts \^a) {i'j 



Ah 



where \i) — I^a^s)- Using Eq. (19), we obtain 

(21) p.«diag(^,^,...,^)+f7,, 

where il,A is a zero diagonal matrix with matrix elements of O [^JNb/N <C 1/Na) (sum 
oi Nb ^ 1 terms of order 1/N with random phases). Neglecting Ua in (21), the reduced 
von Neumann entropy of subsystem A is given by S{pa) ~ log(^A), the maximum 
entropy that the subsystem A can have. 

3'1. Page's formula. - The exact mean value {EabHiP})) of the bipartite entanglement 
is given by Page's formula, obtained [16] by considering the ensemble of random pure 
states drawn according to the Haar measure on U{N): 



(22) Sp = {EABiim = {S{pa)) = {S{pb)) = logNA 



Na 



2Nb In 2 ' 



where ( • ) denotes the (ensemble) average over the uniform Haar measure. For logNA S> 
1, (Eab) is close to its maximum value Sjs^'^d^)) — fog-^A ^ 1- Note that, if we fix 
Na and let Nb ^ oo, then (Eab) tends to E^^^. 

Remarkably, if we consider the thermodynamic limit, that is, we fix Na/Nb and let 
Na -^ oo, then the reduced von Neumann entropy concentrates around its average value 
(22). This is a consequence of the so-called concentration of measure phenomenon: the 
uniform measure on the /c-sphere S'^ in ]R'^+^ (parametrized, for instance, by fc = A^^ — 2 
angles in the Hurwitz parametrization [39, 24]) concentrates very strongly around the 
equator when k is large: Any polar cap smaller than a hemisphere has relative volume 
exponentially small in k. This observation implies, in particular, the concentration of the 
entropy of the reduced density matrix pA around its average value [40, 18]. This in turn 
implies that when the dimension N of the quantum system is large it is meaningful to 
apply statistical methods and discuss typical (entanglement) behavior or random states, 
in the sense that almost all random states behave in essentially the same way. 

3'2. Lubkin's formula. - For random states, the average value of the purity of the 
reduced density matrices pA and pB is giben by Lubkin's formula [41]: 

(23) Pr.^{PipA)) = {PipB)) ^^ + ^^ 



NaNb + 1 
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Nothe that, if we fix Na and let Nb -^ oo, then P^ tends to its minimum value 1/iV^. 
If we fix Na/Nb and let Na -^ oo, then Pl -^ 0. For large N, the variance 



(24) 4 = (P2) _ Pi 



N2-' 



so that the relative standard deviation 



(25) 



V2 



Pl Na + Nb 

tends to zero in the thermodynamic limit Na -^ oo (at fixed Na/Nb)- For balanced 
bipartition, corresponding to Na = Nb = vN, we have 

(26) ^=o(4= 

^ ' Pl \Vn 

Note that the fact that crp / Pl -^ when A^ ^ oo is again a consequence of the concen- 
tration of measure phenomenon [40, 18]. A derivation of Eqs. (23) and (24) is presented 
in Appendix B. 



4. Pseudo-random states 

The generation of a random state \^) is exponentially hard. Indeed, starting from 
a fiducial n^-qubit state |0) = |0. . .0) one needs to implement a typical (random) uni- 
tary operator U (drawn from the Haar measure on U{N = 2"')) to obtain |^) = C/|0). 
Since U is determined by 4"« — 2 real parameters (for instance, the angles of the Hur- 
witz parametrization [39, 24]), its generation requires a sequence of elementary one- and 
two-qubit gates whose length grows exponentially in the number of qubits. Thus, the 
generation of random states is unphysical for a large number of qubits. On the other 
hand, one can consider the generation of pseudo-random states that could reproduce the 
entanglement properties of truly random states [22, 23, 24, 25, 26, 27]. In Refs. [25, 26] 
it has been proven that the average entanglement of a typical state can be reached to a 
fixed accuracy within 0('rvi) elementary quantum gate. This proof holds for a random 
circuit such that U is the product, 

(27) U = WtWt^i---W2Wi, 

of a sequence of t = 0{nT) two-qubit gates Wk independently chosen at each step as 
follows: 

• a pair of integers (i, j), with i ^ j \s chosen uniformly at random from {1, ..., n^} 



9/' 

• singlc-qubit gates (unitary transformations) V^''' and V^'" ^ drawn independently 



/(*) ^^A T/O) 

from the Haar measure on C/(2), are applied; 
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• a CNOT''*'-'' gate with control qubit i and target qubit j is applied. (^) 
Therefore, 

(28) Wk = CNOT^^-'V^^'V^^^^ 

The proof [25, 26] that (27) generates to within any desired accuracy the entanglement 
of a random state in a polynomial number of gates is based on the fact that the evolution 
of the purity of the two subsystems A and B can be described following a Markov chain 
approach, with gap in the Markov chain given by A(ng) > p{nq), where piriq) — 0{n~Z^). 
If the density matrix pt of the overall system is expanded in terms of Pauli matrices. 



(29) Pt= Y. 

ao,---,an„-l=0 



Jao,...,a„,-i) (ao) ^ (a„,_i) 



where 



(30) cr'--""'-^^ = ^Tr(.r^«-«<-v"A 



N' 



with (To ^ -f, cTi = cTj;, (T2 = (T„, (73 == (Tz, wc obtaiu 



(31) P{pA,) = P{PB,) = N^Nl J2 [c|-'- — '°-°y . 

As usual in this paper, we consider ua + ns = nq and NaNb = N . In the case of model 
(28) , the evolution of the column vector 

(32) c,^^*[(c^-°))^(c^-"•^))^...,(c^-^))^] 



is described [25, 26] by a Markov chain: 

(33) c1+^=Mc1. 



{^) By definition, the CNOT*'-^' gate acts on the states {|a;y) = \x)i(Si\y)j = jOO), jOl), jlO), jll)} 
of the two-qubit computational basis as follows: CNOT*'-^' turns jOO) into jOO), jOl) into jOl), 
jlO) into jll), and jll) into jlO). The CNOT (controUed-NOT) gate flips the state of the 
second (target) qubit if the first (control) qubit is in the state jl) and does nothing if the first 
qubit is in the state |0). In short, CNOT(''J'|x)i|j/)j = \x)i\x®y)j, with x,y £ {0,1} and 
© indicating addition modulo 2. By definition, given input bits x, y, the XOR gate outputs 
i (B j- Therefore, the (quantum) CNOT gate acts on the states of the computational basis 
as the (classical) XOR gate. However, the CNOT gate, in contrast to the XOR gate, can 
also be applied to any superposition of the computational basis states. The CNOT gate is 
the prototypical two-qubit gate that is able to generate entanglement. For instance, CNOT 
maps the separable state |i/>) = -TsdO) + |1))|0) onto the maximally entangled (Bell) state 

ji^"*") = CNOT 1 1/)) = -i=(|00) + jll)). Any unitary operation in the Hilbert space of n, qubits 
can be decomposed into (elementary) one-qubit and two-qubit CNOT gates. 
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Therefore, the asymptotic decay of purity is determined by the second largest eigenvalue 
1 — A(ng) of the matrix M (the largest eigenvalue is the unit eigenvalue): 

(34) \{P{pA.t)) ~ Pl\ -{I- Ain,)Y = exp{[ln(l - A(n,))]i}. 

For model (28), one obtains [25, 26] A{nq) > p{nq), with p{ng) = 0{n~'^). 

Alternatively, a two-qubit gate different from CNOT [27] or Wt chosen from the U{A) 
Haar measure and acting on a a pair i,j of qubits (i ^ j) randomly chosen at each step, 
can be used [27, 28]. Numerical results for these models [25, 26, 27] indicate that the 
above analytic bound is not optimal: there is numerical evidence that A(nq) = 0(n~^) 
and therefore 0{ni) steps are sufficient to generate the average entanglement of a random 
state. A different strategy to efficiently approximate the entanglement content of random 
states is based on quantum chaotic maps [24, 42] and will be discussed in Sec. 6. 

Finally, we note that the Markov chain approach allows an easy derivation of Lubkin's 
formula (23). Let us consider Wk chosen from the U{4) Haar measure. In this case, the 
Markov matrix 

with M}- acting non trivially (differently from identity) only on the subspace spanned by 
qubits i and j {i, j — I, ..., Hq and i ^ j). After averaging over the uniform Haar measure 

on t/(4) one can see [27] that Mj- preserves identity (cr^ (g) crj — > o-^ (g) crj ) and 

uniformly mixes the other 15 products <7^ (g) ay' [27]. Matrix elements are therefore 

[M^f]o,o = 1, [M^j\,x = [M^j^Uo = 0, and [M^'^],,,, = 1/15 for x,x' e {1, ..., 15}. 

The matrix M has an eigenvalue equal to 1 (with multiplicity 2) and all the other 
eigenvalues smaller than 1. The eigenspace corresponding to the unit eigenvalue of matrix 
M is spanned by the column vectors 

(36) «o-*(l,0,...,0), i;i=*(0,l,...,l). 

The asymptotic equilibrium state c^(oo) = limt^oo c^(0 i'' however uniquely determined 
by the constraints Tr(pt) — 1 and Tr(p() — 1, which impose 

(37) .0 - [c^-y - ]^' -^ - E [cl—- -)]^ = ^. 

ao,...,ar.,-l#(0,...,0) 

Finally, we obtain 

2/ N 1 If/ N- 1 N- 1 , 

(38) c (oo) = xovo + a^i^^^— Y«i = ^ (1, ]^^— Y' -' ]7^— y) 

and, after substitution of the components of this vector into Eq. (31), 



2 



(39) PiPA.) = ^ 



1 + (^1-1)^ 



which immediately leads to Lubkin's formula (23). 
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5. Multipartite entanglement 

The characterization and quantification of multipartite entanglement is a challenging 
open problem in quantum information science and many different measures have been 
proposed [2, 3]. To grasp the difficulty of the problem, let us suppose to have n parties 
composing the system we wish to analyze. In order to obtain a complete characterization 
of multipartite entanglement, we should take into account all possible non-local corre- 
lations among all parties. It is therefore clear that the number of measures needed to 
fully quantify multipartite entanglement grows exponentially with the number of qubits. 
Therefore, in Ref. [43] it has been proposed to characterize multipartite entanglement 
by means of a function rather than with a single measure. The idea is to look at the 
probability density function of bipartite entanglement between all possible bipartitions 
of the system. For pure states the bipartite entanglement is the von Neumann entropy of 
the reduced density matrix of one of the two subsystems: Eab{\iP)) — S{pa) = S{pb)- 

It is instructive to consider the smallest non-trivial instance where multipartite en- 
tanglement can arise: the three-qubit case. Here we have three possible bipartitions, 
with ua = I qubit and ub = 2 qubits. For a GHZ state [44], 

(40) |GHZ) = i=(|000) + |lll)), 

we obtain pa ~ 2 ^or all bipartitions, and therefore 

(41) p{Eab) = Se^bA, 

namely there is maximum multipartite entanglement, fully distributed among the three 
qubits. (^) Note that in this case pB — |(|00)(00| + |11)(11|) is separable and therefore 
the pairwise entanglement between any two qubits is equal to zero. 
For a W state [46] , 

(42) |W)--^(|100) + 1010) + 1001)), 

we obtain pA = ||0)(0| + i|l)(l| for all bipartitions, and therefore 

(43) p{Eab) = Se^^^b, 

where E = —^ log (|) — i log (i) « 0.92, namely the distribution p{Eab) is peaked but 
the amount of multipartite entanglement is not maximal. 
As a last three-qubit example, let us consider the state 

(44) 1^) = -^(jOOO) + 1110)) = i=(|00) + |11))|0), 



(*) The problem of finding, for a generic number n, of qubits, maximally multipartite entan- 
gled states, that is, pure states for which the entanglement is maximal for each bipartition, is 
discussed in [45]. 
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where the first two qubits are in a maximally entangled (Bell) state, while the third one is 
factorized. In this case, pA = | if subsystem A is one of the first two qubits, pA = |0)(0| 
otherwise. Hence, 

2 1 

(45) p(Eab) = g^B^B.i + ^Seab,o, 

namely the entanglement can be large but the variance of the distribution p{Eab) is also 
large. 

For sufficiently large systems {N — 2"« 3> 1), it is reasonable to consider only bal- 
anced bipartitions, i.e., with ua = ns {riA + ns — riq), since the statistical weight of 
unbalanced ones {ua ^ ns) becomes negligible [43]. If the probability density has a large 
mean value (Eab) ^ nq (( • ) denotes the average over balanced bipartitions) and small 
relative standard deviation (Jab/{Eab) ^ 1, we can conclude that genuine multipartite 
entanglement is almost maximal (note that Eab is bounded within the interval [0,nq]). 
As I shall discuss in Sec. 6, this is the case for random states [43] (see also Ref. [47]).(®) 

6. Quantum chaos map and entanglement 

The relation between chaos and entanglement is discussed, for instance, in [49, 50, 51, 
52, 53]. More specifically, the use of quantum chaos for efficient and robust generation 
of pseudo-random states carrying large multipartite entanglement is nicely illustrated by 
the example of the quantum sawtooth map [54, 36, 7]. This map is described by the 
unitary Floquet operator U: 

(46) IV't+i) = U\A) = e-^"V2 ^.fc(«-.)V2|^^^^ 

where h = —id/dO, [9^fi] = i (wc set h ^ I) and the discrete time t measures the 
number of map iterations. In the following I will always consider map (46) on the torus 
< ^ < 27r, —TT < p < TT, where p = Tn. With an rig-qubit quantum computer we are 
able to simulate the quantum sawtooth map with N = 2"'' levels; as a consequence, 9 
takes N equidistant values in the interval < 9 < 2tt, while n ranges from —N/2 to 
N/2 — 1 (thus setting T = 2Tr/N). Wc are in the quantum chaos regime for map (46) 
when K = kT > ov K < —4; in particular, in the following I will focus on the case 
K ^ 1.5. 

There exists an efficient quantum algorithm for simulating the quantum sawtooth 
map [54, 7]. The crucial observation is that the operator U in Eq. 46 can be written as 

the product of two operators: JJk — e*'^'^^^^' /^ and Ut — e^*"^" /^, that are diagonal 
in the 9- and in the n-representation, respectively. Therefore, the most convenient way 
to classically simulate the map is based on the forward-backward fast Fourier transform 
between 9 and n representations, and requires O(A^logA^) operations per map iteration. 
On the other hand, quantum computation exploits its capacity of vastly parallelize the 
Fourier transform, thus requiring only 0((logA^)^) one- and two-qubit gates to accom- 
plish the same task [54, 7]. In brief, the resources required by the quantum computer to 



(^) For a Tiq-qubit GHZ state, |GHZ) = -i=(|0...0) -I- |1. ..!)), the distribution p(Eab) is peaked 
but at a small and essentially riq-independent value, while for the cluster states [48] the average 
entanglement {Eab) is large and increases with Uq but also the variance is large [43]. 
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Fig. 1. - Time evolution of the average bipartite entanglement of a quantum state, starting from 
a state of the computational basis (eigenstate of the momentum operator n), and recursively ap- 
plying the quantum sawtooth map (46) at K = 1.5 and, from bottom to top, n, = 4, 6, 8, 10, 12. 
Dashed lines show the theoretical values of Eq. (47). Inset: convergence of {EAB){t) to the 
asymptotic value {E^g"^) in Eq. (47); time axis is rescaled with l/n,. This figure is taken from 
Ref. [421. 



simulate the sawtooth map are only logarithmic in the system size N^ thus admitting an 
exponential speedup, as compared to any known classical computation. The sawtooth 
map and the quantum algorithm for its simulation are discussed in details in Appendix C. 
Let us first compute the average bipartite entanglement (Eab) as a function of the 
number t of iterations of map (46). Numerical data in Fig. 1 exhibit a fast convergence, 
within a few kicks, of this quantity to the value 



(47) 



(E. 



AB 



Uq 

2 



1 



21n2 



expected for a random state according to Page's formula [16] (note that this result is 
obtained from Eq. (22) in the special case n^ — n^). Precisely, as shown in the inset of 
Fig. 1, {Eab) converges exponentially fast to {E^^'^), with the time scale for convergence 
oc Uq. Therefore, the average entanglement content of a true random state is reached 
to a fixed accuracy within 0{ng) map iterations, namely 0{n^) quantum gates. I stress 
that in our case a deterministic map, instead of random one- and two-qubit gates as in 
Ref. [25, 26, 27, 28], is implemented. Of course, since the overall Hilbert space is finite, 
the above exponential decay in a deterministic map is possible only up to a finite time 
and the maximal accuracy drops exponentially with the number of qubits. I also note 
that, due to the quantum chaos regime, properties of the generated pseudo-random state 
do not depend on initial conditions, whose characteristics may even be very far from 
randomness (e.g., simulations of Fig. 1, start from a completely disentangled state). 

As discussed above, multipartite entanglement should generally be described in terms 
of a function, rather than by a single number. I therefore show in Fig 2 the probabil- 
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Fig. 2. - Probability density function of the bipartite von Neumann entropy over all balanced 
bipartitions for the state \tpt), after 30 iterations of map (46) at K = 1.5. Various histograms 
are for different numbers of qubits: from left to right Uq = 8, 10, 12; dashed curves show the 
corresponding probabilities for random states. Inset: relative standard deviation oab/ {Eab) as 
a function of n, (full circles) and best exponential fit (tab /{Eab) ~ g-O-^Sn, (continuous line); 
data and best exponential fit <jab/{Eab) ~ e~"'" for random states are also shown (empty 
triangles, dashed line). This figure is taken from Ref. [42]. 



ity density function p{Eab) for the entanglement of all possible balanced bipartitions 
of the state |f/'t=3o)- This function is sharply peaked around {E^^^), with a relative 
standard deviation uab/ (Eab) that drops exponentially with n, (see the inset of Fig. 2) 
and is small (~ 0.1) already at riq = 4. For this reason, we can conclude that multi- 
partite entanglement is large and that it is reasonable to use the first moment (Eab) 
oi p{Eab) for its characterization. The corresponding probability densities for random 
states is also calculated (dashed curves in Fig. 2); their average values and variances are 
in agreement with the values obtained from states generated by the sawtooth map. As 
we have remarked in Sec. 3, the fact that for random states the distribution p{Eab) 
is peaked around a mean value close to the maximum achievable value E'^^ = nq/2 
is a manifestation of the concentration of measure phenomenon in a multi-dimensional 
Hilbert space [40, 18]. 



7. Stability of multipartite entanglement 

In order to assess the physical significance of the generated multipartite entanglement, 
it is crucial to study its stability when realistic noise is taken into account. Hereafter 
I model quantum noise by means of unitary noisy gates, that result from an imperfect 
control of the quantum computer hardware [55]. The noise model of Ref. [56] is fol- 
lowed. One-qubit gates can be seen as rotations of the Bloch sphere about some fixed 
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axis; I assume that unitary errors slightly tilt the direction of this axis by a random 
amount. Two-qubit controlled phase-shift gates are diagonal in the computational basis; 
I consider unitary perturbations by adding random small extra phases on all the compu- 
tational basis states. Hereafter I assume that each noise parameter e^ is randomly and 
uniformly distributed in the interval [—£,+&]; errors affecting different quantum gates 
are also supposed to be completely uncorrelated: every time we apply a noisy gate, noise 
parameters randomly fluctuate in the (fixed) interval [— e, -|-£]. 

Starting from a given initial state |V'o)j the quantum algorithm for simulating the 
sawtooth map in presence of unitary noise gives an output state |V'e/,t) that differs from 
the ideal output JV't)- Here ej = {ei,e2, ■■■,£nd) stands for all the Ud noise parameters 
Ei, that vary upon the specific noise configuration (rid is proportional to the number of 
gates). Since we do not have any a priori knowledge of the particular values taken by the 
parameters Si, the expectation value of any observable A for our n^-qubit system will be 
given by Tr[pg^tA], where the density matrix p^^t is obtained after averaging over noise: 

(48) P''={y)j ' j dei\i^e,,t){^e,A- 

The integration over sj is estimated numerically by summing over M random realizations 
of noise, with a statistical error vanishing in the limit M -^ oo. The mixed state p^ may 
also arise as a consequence of non-unitary noise; in this case Eq. (48) can also be seen 
as an unraveling of p^ into stochastically evolving pure states |V'e/)j each evolution being 
known as a quantum trajectory [57, 58, 59, 60]. 

I now focus on the entanglement content of pe^t- Unfortunately, for a generic mixed 
state of Uq qubits, a quantitative characterization of entanglement is not known, neither 
unambiguous [2, 3]. Anyway, it is possible to give numerically accessible lower and upper 
bounds for the bipartite distillahle entanglement Ej^J(p^): 

(49) max{5(pe,A) - Sip,),0} < Efj^ipe) < log||p^«|| , 



where pe^A = T^'csipe) and \\p'^^\\ = TrJ {pe^)'^ Pe^ denotes the trace norm of the partial 
transpose of pe with respect to party B (see Appendix A'l for the definition of the partial 
transposition operation) . 

In practice, the quantum algorithm for the quantum sawtooth map is simulated in 
the chaotic regime with noisy gates and the two bounds in Eq. (49) for the bipartite 
distillable entanglement of the mixed state p^^ti obtained after averaging over M noise 
realizations, are evaluated. A satisfactory convergence for the lower and the upper bound 
is obtained after M ^ \/N and J\f ^ N noise realizations, respectively. The first moment 
of the lower {Em) and the upper [Em] bound for the bipartite distillable entanglement 
is shown as a function of the imperfection strength in Fig. 3, upper panels. The various 
curves are for different numbers Uq of qubits; M depends on rig and is large enough to 
obtain negligible statistical errors (smaller than the size of the symbols). The relative 
standard deviation of the probability density function (over all balanced bipartitions) 
for the bipartite distillable entanglement is shown in the lower panels of Fig. 3. Like 
for pure states, we notice an exponential drop with n^; the distribution width slightly 
broadens when increasing imperfection strength e. We can therefore conclude that an 
average value of the bipartite bipartite distillable entanglement close to the ideal case 
e = implies that multipartite entanglement is stable. 
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Fig. 3. - Upper graphs: lower {Em) (left panel) and upper bound {Em) (right panel) for the 
bipartite distillable entanglement as a function of the noise strength at time t = 30. Various 
curves stand for different numbers of qubits: Uq = A (circles), 6 (squares), 8 (diamonds), 10 
(triangles up), and 12 (triangles down). Lower graphs: relative standard deviation of the prob- 
ability density function for distillable entanglement over all balanced bipartitions as a function 
of Uq, for different noise strengths e. Dashed lines show a behavior a/ {E) ~ e""''^ and are 
plotted as guidelines. This figure is taken from Ref. [42]. 



In order to quantify the robustness of multipartite entanglement with the system size, 
let us define a perturbation strength threshold e^^' at which the distillable entanglement 
bounds drop by a given fraction, for instance to 1/2, of their e = value, and analyze 
the behavior of e'^^^ as a function of the number of qubits. Numerical results are plotted 
in Fig. 4; both for lower and upper bounds we obtain a power-law scaling close to 



(50) 



AR) 



l/riq 



It is possible to give a semi-analytical proof of the scaling (50) for the lower bound 
measure, that is based on the quantum Fano inequality [61], which relates the entropy 
S{pe) to the fidelity F = {iPt\peAi't)- 



(51) 



S{Pe) < HF) + (1 - F) \og{N' - 1), 



where h{x) 



F 



-JS Ugt 



— a::log(a::) — (1 
[56, 62], with 7 



x) log(l — x) is the binary Shannon entropy. Since 
■• 0.28 and Ug — 3n^ + Ug being the number of gates 



18 



G. BENENTI 





. 1 


1 ' 1 ' 1 ' 1 . 


0.010 




- 


(R) 




- 




■ 




m 




•^ 


0.008 


- \ 


^s^ 


0.007 


- 


^* 






^x ^^ 


0.006 


- 


\ 






■^^ 's 


0.005 


- 


\" 






\ 


0.004 


" 


^^> 




1 


1 ,1.1,1 



10 12 



0.028 


"•^ 


' 1 ' 


1 ' 1 ' 


_ 


(R) 


X 


• ;=15 


- 


S/f 


■ 


^^ 


f = 30 


- 


M 


\^ 




- 


0.020 




• 




- 


'^.^ 


0.016 


- 


"•s 




- 




0.012 










>^ 


0.008 


- 


1 ,1,1,1 



10 12 



Fig. 4. - Perturbation strength at which the bounds of multipartite entanglement halve (lower 
bound on the left panel, upper bound on the right panel), as a function of the number of 
qubits. Dashed lines are best power-law fits of numerical data: e 



(R) 



at t — 15, 



,(fl) 



at f = 30, for both lower and upper bounds. This figure is taken from Ref. [42]. 



required for each map step, we obtain, for e^Ugt <C 1, 



(52) 



S{pe) < le^ngt 



log(7e^%t) + 2"? + j^ 



For sufficiently large systems the second term dominates (for Uq — 12 qubits, i = 30 and 
e ~ 5 X 10~^ the other terms are suppressed by a factor ~ 1/10) and, to a first approx- 
imation, we can only retain it. On the other hand, an estimate of the reduced entropy 
S{pa) is given by the bipartite entropy (47) of a pure random state [16]. Therefore, 
from Eq. (49) we obtain the following expression for the lower bound of the distillable 
entanglement: 



(53) 



<B^(Pe) > ^ 



1 



21n2 



Q-^TVie^t . 



{D), 



From the threshold definition E)^^ (p, 
that is valid when n„^ 1: 



AD) 



k^AsiPo) = ^S{pa) we get the scaling (50), 



(54) 



4^) ^ 1/^247^2 t. 



Notice that, for small systems as the ones that can be numerically simulated (see data in 
Fig. 4), the first term of Eq. (52) may introduce remarkable logarithmic deviations from 
the asymptotic power-law behavior. At any rate, the scaling derived from Eq. (53) is in 
good agreement with the above shown numerical data, and also reproduces the prefactor 
in front of the power-law decay (50) up to a factor of two. 
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8. Detecting entanglement of random states 

The entanglement content of high-dimensional random pure states is almost maximal: 
nevertheless, in this section I will demonstrate that, due to the complexity of such states, 
the detection of their entanglement is rather difficult [63] . 

8'1. Random states and the quantum to classical transition. - A motivation to the 
study of the detection of random states comes from considerations of the quantum to 
classical transition. Since random states carry a lot of entanglement and entanglement 
has no analogue in classical mechanics, one can immediately conclude that random states 
are highly non-classical. On the other hand, as we have seen in Sec. 6, pseudo-random 
states with properties close to those of true random states can be efficiently generated 
by dynamical systems (maps) in the regime of quantum chaos. In such chaotic maps 
the classical limit is recovered when the number of levels N ^ oo. Therefore, one can 
argue that for large random states, i.e., in the limit A^ — > oo, the quantum expectation 
value of an operator with a well defined classical limit will be close to its classical micro- 
canonical average. According to this picture random states in a way "mimic" classical 
microcanonical density. Expectation values are therefore close to the classical ones. 

At first sight this is in striking contrast with the almost maximal entanglement of 
such states. However, as I shall discuss in the following, the contradiction is only ap- 
parent [63] . The detection of entanglement for a random state appears very difficult at 
large A^, as it would demand the control of very finely interwoven degrees of freedom 
and a measurement resolution inversely proportional to N , which seems hardly feasible 
experimentally. Therefore, as far as the detection of entanglement is concerned, high 
dimensional random states are effectively classical. {^'^) 

Moreover, coarse graining naturally appears. For instance, one could repeat several 
times the measurement of an entanglement witness (the definition of entanglement wit- 
ness is provided in Appendix A'2) for a random state and the prepared random state 
would be different from time to time due to unavoidable experimental imperfections. Let 
us model this problem by considering mixtures of m pure random states, namely 



(55) p = Y,-\^^){^, 



m 

where the \^i) are mutually independent random pure states, but in general they are not 
orthogonal. I am going to show that the detection of entanglement is even more difficult 
for these mixed states, as it requires a number of measurements growing exponentially 
with TO. 

It is interesting to remark that there are other physical contexts in which formally 
the same kind of coarse graining naturally appears: 

• (i) Time averaging - For example, if a state [■(/') undergoes a time evolution |V'(^)) — 
U{t)\tl?j given in terms of some unitary dynamics U{t), then the time average of a 
physical observable A over an interval T is given by the expectation value Tv{Ap) 



(^'^) Of course, this remark does not call into question the utility of high dimensional random 
states for quantum information processing. 
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for the mixed state 

(56) p=^J dt|^(t))(^(i)|, 

which has an effective rank m « T/icorrj where icorr is a dynamical correlation time 
of the dynamics U{t). For a quantum chaotic evolution U{t), the state |'0(i)) can 
be, after some time, arguably well described by a random state and the correlation 
time icorr is expected to be short, so p in Eq. (56) may be well approximated by a 
mixture of m uncorrelated random states analogous to Eq. (55). 

• (ii) Phase space averaging - Sometimes it is useful to represent quantum states in 
terms of distribution functions in the classical phase space, like the Husimi function 
(see, e.g., [31]), which can be understood as a convolution of the Wigner function 
or its coarse graining over a phase space volume 27rft (to simplify writing, let us 
consider systems with one degree of freedom). In fact, the Husimi function of a 
pure state can be understood as a Wigner function of the following mixed state: 



(57) p= -^^ I dqdp cxp 



— (ag^ +a ^/) 



T{q,pM{i;\T\q,p), 



where T(q,p) are unitary phase space translation operators, and a is an arbitrary 
(squeezing) parameter. A random pure state {ip) has a Wigner function with ran- 
dom sub-Planck structures with phase space correlation length Icon ^ ft which is 
semi-classically smaller than the coarse-graining width ^ h^^^, so p in Eq. (57) 
can be again considered as a mixture (Eq. (55)) of m random pure states with 

8'2. Unknown random states. - In this section is is assumed that the random state \^) 
whose entanglement we would like to detect is unknown so that we are not able to use an 
optimal entanglement witness W for a particular \ip). The best one can do is to choose 
some fixed witness W in advance, independently of the state. Since I am interested in the 
average behavior over unitary invariant ensemble of pure random states, W can be chosen 
to be random as well. That is, in the present section I am going to study detection of 
entanglement with a random entanglement witness, whose precise definition will be given 
later. What I want to calculate is the distribution of the expectation values (V'lVFlV') for 
a fixed W and an ensemble of random pure states \ip). Averaging over random states jf/') 
we see that the average expectation value (V'lM^IV') is 

(58) f dr{'4j\W\ijj)^Tr{W)/N, 



where • = J dV* denotes an integration over the C/(7V)-invariant (Haar) distribution 
of pure states |^), and I used the fact that for a random state |V') = X^i Ci|i) we have 
CiC* — Sij /N. Let us fix normalization of the entanglement witness W such that Tr(VF) — 

1. Therefore, the average expectation value (■(/'IVFI^) scales ex 1/A^. It is therefore 
convenient to define the rescaled quantity w = N{ip\W\'ip) such that uJ = 1, independently 
of the dimension N. 
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Here and in the following, the investigation is limited to decomposable entanglement 
witnesses (see Appendix A'2) of the form W = Q"^^ , with Q positive semidefinite oper- 
ator. 

I first consider the case when Q is a simple rank one projector, that is, W is given 
hy W — {\<l)) {(j)\)^^ . If \(f>) is a state with a large Schmidt number r ^ VN, as it is 
typical for random |(/)) (I consider equal size subsystems, Na — Nb — Vn), then one 
can show [63] that the probability density p{w) = dV/dw converges to a Gaussian in the 
limit N ^ oo, 

(59) p{w) = -^cxp{-{w-lf/2). 



Numerical results for finite A^ = 2^" are shown in Fig. 5 (top). The probability of 
measuring negative w, i.e., of detecting entanglement, is J_ p{w)dw and therefore 



(60) V{w < 0) = (1 - erf(l/V2))/2 w 0.159. 

Note that this entanglement detection probability is independent of the details of \(p), 
provided that its Schmidt number r is large, more precisely r (x N. 

Since it appears difficult to measures witness operators corresponding to states |0) 
with large Schmidt number r, it is interesting to consider the opposite limit of small r. 
In particular, let us consider the extreme case of rank r = 2. We therefore have only 
two nonzero terms in the Schmidt decomposition (10), corresponding to the nonzero 
eigenvalues, pi = X and p2 = 1 ^ ^^j of the reduced density matrix pA = Ti b {\4') {(f>\) ■ In 
this case, we obtain [63] 

(a^{Ae-f + (l-A)e-T-x} + __i=_e-7^S=^ : w > 0, 



(61) p{w) 



1 



4^A(l-A)+2 



ev^<i^^ : w <0. 



Results of numerical simulation for p{w) for two cases, A = 1/2 and A = 1/26, are 
compared in in Fig. 5 (bottom). The probability of detecting entanglement, i.e., of 
measuring negative values of w is 

(62) r{w < 0) = 1/(4 + 2/v/A(l - A)). 

For instance, V{w < 0) = 1/8 = 0.125 when A ^ 1/2 and r{w < 0) = 5/72 w 0.07 when 
A ~ 1/26. Note that in the limit A —> 0, i.e., of a pure separable state for \(p), w is always 
positive with an exponential distribution. 

I emphasize that, neglecting problems related to finite measurement resolution (an 
issue that will be discussed in Sec. 8'3), the entanglement detection probability is, for 
pure states, independent of the system size A^. 

So far I have discussed only the case when Q is of rank one, Q = \4>){4'\- In general, 
one can consider Q of rank k: Q — ^^ di\(f>i){(pi\. Since I have fixed Tr(VF) = 1, then 
^j di = 1. Assuming for simplicity that all di are the same, di — 1/fc, we obtain [63], in 
the limit N —^ oo, 

(63) p{w) ^ yAe-'=(-i)V2. 
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Fig. 5. - Distribution of w = Ai'{i/>| Wl-i/;) for random vectors \ip) and a single 14^ = (|0)(0j) ^ 
with a random |<jf)) (top) and with |0) having two nonzero Schmidt coefficients vA and \/l — A 
(two cases are shown, A = f/2 and A = 1/26) (bottom). Histograms correspond to numerical 
simulations for A'' = 2^", curves to the theoretical predictions of Eqs. (59) (top) and and (61) 
(bottom). This figure is taken from Ref. [63]. 



Therefore, the entanglement detection for pure states is more efficient if one consider 
W = Q^'^ , with Q rank-one projector. 

It is also interesting to consider the case in which Q is of rank 1 but we wish to detect 
the entanglement of mixed states, for instance of a mixture of m pure random states, as 
given in Eq. (55). In this case, the resulting distribution p{w) is, in the limit A^ -^ oo, a 
Gaussian of variance 1/m [63]. Because the distribution p{w) becomes narrowly peaked 
about its mean W = 1 with increasing m, the probability of measuring negative values 
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decreases with m, that is, we obtain 



(64) n. < 0) = i^i^^!^ , vt^""^^- 

This probabihty decays to zero exponentially with m. Therefore, the detection of entan- 
glement for a mixture of random states is very hard. This result outlines the importance 
of coarse graining to explain the emergence of classicality. For random pure states, a fi- 
nite success probability in the detection of entanglement exists also in the limit in which 
the Hilbert space dimension N ^ oo. This implies that chaotic dynamics alone is not 
sufficient to erase any trace of entanglement when going to the classical limit, provided 
that ideal measurements are possible. On the other hand such erasure becomes very 
efficient when coarse graining is taken into account, for instance when mixtures instead 
of pure states are considered. 

8'3. Known random states. - In this section it is assumed that the random state |V') 
whose entanglement we want to measure is known in advance and furthermore, that we 
are able to prepare an arbitrary decomposable entanglement witness. In addition, we have 
to assume that our state \ip) is neither separable, nor bound entangled (see Appendix A), 
which is true with probability which converges to one exponentially in N. Therefore, 
for each \ip) we can prepare an optimal entanglement witness, such that its expectation 
value will be minimal. As far as decomposable entanglement witnesses are concerned, the 
optimal choice of H^ = VFopt is to take for Q a projector to the eigenspace corresponding 
to the minimal (negative) eigenvalue Amin oi p^^ , Wopt = {\(f>min) {4'min\)'^^ ■ The maximal 
violation of positivity is therefore 

(65) Tr(Woptp) = -|A„,i„(pT^)|. 

If we are able to measure the entanglement witness Wopt with a given precision it is the 
size of Amin which determines the difficulty of detecting entanglement in | ■;/')■ Note that 
the optimal entanglement witness Wopt depends on the state \tp). For each state |V') we 
have to pick a different Wopt • 

The expectation value of the minimal eigenvalue equals Amin = —4:/yN [63]. In 
fact, the distribution of Amin becomes strongly peaked around — 4/viV with diminishing 
relative fluctuations as A'^ -^ oo. When we mix several independent (in general non- 
orthogonal) random vectors, p — ^j \4'i){'4'i\/'n^i the minimal eigenvalue Amin increases 
and the distribution becomes increasingly sharply peaked (for m —^ oo we get p —^ 1/N 
with all eigenvalues being equal to 1/N). Note that the average minimal eigenvalue Amin 
is positive for m > to*, with to* « AN. 

Although von Neumann entropy of a random state is large all eigenvalues of p^^ are 
very small and will therefore be hard to detect. If we assume that we are able to measure 
values of |Tr(W/9)| < e then we can, depending on the scaling of e with N, tell for which 
values of to, the detection of entanglement is possible. If e does not depend on N, i.e., 
precision does not increase with N, then for sufficiently large N, such that A/^/N < e, 
detection of entanglement will be impossible. Already a single random state becomes 
from the viewpoint of entanglement detection "classical", since measuring a negative 
expectation value of its optimal entanglement witness is below the detection limit. If on 
the other hand we are able to measure e which decreases as 1/yN, the critical TO-crit, 
beyond which the entanglement detection is impossible, will be independent of N , i.e.. 
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in the limit iV — s- oo the ratio merit/ \ N ~^ 0. If however we are able to detect very 
small expectation values of order 1/A'^, then Wcrit will be proportional to N . Furthermore, 
even with arbitrary accuracy, detection of entanglement with decomposable entanglement 
witnesses is in practice impossible beyond m = m* ex N (that is, when Amin becomes 
positive) . 

I would like to stress once more than the detection difficulties are a consequence 
of the complexity of random states. If instead one considers "regular states" such as 
the GHZ state [44] |GHZ) = -ij(|0...0) + |l...l)), then the optimal witness is Wopt = 
(|'/'mi„)('/'mi„|)'^^, with |,?;.„iin) - ^ (|0.. .0) 1 1.. .1) - 1 1.. .1) |0.. .0) ) which corrcspouds to the 

minimal eigenvalue Amin = ^1/2 of (|GHZ)(GIIZ|)'^^. Since the value of Amin is —1/2 
instead of —A/yN as typical for a random state, it turns out that it will be much easier 
to detect entanglement in a "regular" rather than in a random state. This happens in 
spite of the fact that the entanglement content is much larger in a random than in such 
a regular state. 



9. Chaotic environments 

Real physical systems are never isolated and the coupling of the system to the envi- 
ronment leads to decoherence. This process can be understood as the loss of quantum 
information, initially present in the state of the system, when non-classical correlations 
(entanglement) establish between the system and the environment. On the other hand, 
when tracing over the environmental degrees of freedom, we expect that the entanglement 
between internal degrees of freedom of the system is reduced or even destroyed. Decoher- 
ence theory has a fundamental interest, since it provides explanations of the emergence 
of classicality in a world governed by the laws of quantum mechanics [64]. Moreover, it is 
a threat to the actual implementation of any quantum computation and communication 
protocol [7, 8]. Indeed, decoherence invalidates the quantum superposition principle, 
which is at the heart of the power of quantum algorithms. A deeper understanding of 
the decoherence phenomenon is essential to develop quantum technologies. 

The environment is usually described as a many-body quantum system. The best- 
known model is the Caldeira-Leggett model [65, 66, 67], in which the environment is a 
bosonic bath consisting of infinitely many harmonic oscillators at thermal equilibrium. 
More recently, first studies of the role played by chaotic dynamics [68, 69, 70, 71, 72, 73] 
or random environments [74, 75] in the decoherence process have been carried out. 

In the following, it is shown that the many-body environment may be substituted with 
a closed deterministic system with a small number of degrees of freedom, but chaotic [72] . 
In other words, the complexity of the environment arise not from being many-body 
but from having chaotic dynamics. I consider two qubits coupled to a single particle, 
fully deterministic, conservative chaotic "environment" , described by the kicked rotator 
model. It is shown that, due to the system-environment interaction, the entropy of the 
system increases. At the same time, the entanglement between the two qubits decays, 
thus illustrating the loss of quantum coherence. The evolution in time of the two-qubit 
entanglement is in good agreement with the evolution obtained in a pure dephasing 
stochastic model. Since this pure dephasing decoherence mechanism can be derived in 
the framework of the Caldeira-Leggett model [76] , a direct link between the effects of a 
many-body environment and of a chaotic single-particle environment is established. 

Let us consider two qubits coupled to a quantum kicked rotator. The overall system 
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is governed by the Hamiltonian 

(66) H = //(I) + ff (2) + ^(kr) _^ ^(int) ^ 

where H*^*' = uji ax (* = l, 2) describes the free evolution of the two qubits, 

kcos{e)"^d{T-jT) 

3 



the interaction between the qubits and the kicked rotator; as usual, CTq {a = x, y, z) 
denote the Pauli operators for the i-th qubit. Both the cosine potential in ijC^'') and the 
interaction _ff('"*) are switched on and off instantaneously (kicks) at regular time intervals 
T . Let us consider the two qubits as an open quantum system and the kicked rotator 
as their common environment. Note that I chose non- interacting qubits as I want their 
entanglement to be affected exclusively by the coupling to the environment (^^) 

The kicked rotator, as the sawtooth map described in detail in Appendix C, belongs to 
the class of periodically driven systems of Eq. (C.f ), with the external driving described 
by the potential V{9) ~ kcos0, switched on and off instantaneously at time intervals 
T. The evolution from time tT~ (prior to the i-th kick) to time {t + 1)T" (prior to the 
{t + l)-th kick) of the kicked rotator in the classical limit is described by the Chirikov 
standard map: 

, . ( fit+i =nt + A: sin 6*4, 

^ ^ 1 9t+i^et+Tnt+i, 

where (n, 0) are conjugated momentum-angle variables and t = t/T denotes the discrete 
time, measured in number of kicks. As for the sawtooth map, by rescaling n ^ p ~ Tn, 
the dynamics of Eq. (69) is seen to depend only on the parameter K = kT . For K = Q the 
motion is integrable; when K increases, a transition to chaos of the Kolmogorov-Arnold- 
Moser (KAM) type is observed [79, 80]. (^^) The last invariant KAM torus is broken for 
K K, 0.97. li K ^ 1 the phase space is mixed (simultaneous presence of integrable and 
chaotic components). If K increases further, the stability islands progressively reduce 
their size; for K ^ 1 they are not visible any more. In what follows, I always consider 



(^^)As discussed in [70], the chaotic environment model discussed in this section could be 
implemented, at least in principle, using cold atoms in a pulsed optical lattice created by laser 
fields [77] or superconducting nanocircuits [78]. 

(^^) The property of complete integrability is very delicate and atypical as it is, in general, 
destroyed by an arbitrarily weak perturbation that converts a completely integrable system 
into a KAM-integrable system. The structure of KAM motion is very intricate: the motion 
is confined to invariant tori for most initial conditions yet a single, connected, chaotic motion 
component (for more than two degrees of freedom) of exponentially small measure (with respect 
to the perturbation) arises, which is nevertheless everywhere dense. 
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map (69) on the torus < 9 < 27r, — vr < p < n. In this case, the Chirikov standard 
map describes the stroboscopic dynamics of a conservative dynamical system with two 
degrees of freedom which, in the fully chaotic regime K ^ I, relaxes, apart from quantum 
fluctuations, to the uniform distribution on the torus. 
The Hilbert space of the global system is given by 

(70) W = W(i) ® 7^(2) ^ -^(fer) ^ 

where TC'^^^ and Ti^^' are the two-dimensional Hilbert spaces associated to the two qubits, 
and T^C^'") is the Hilbert space for the kicked rotator with TV quantum levels. 

The time evolution generated by Hamiltonian (66) in one kick is described by the 
operator 



(71) 



U^ exp [ - i(fc + ei&P + (7?^)) cos(0')] 

X exp [ — *T^] exp ( — i (5i (Ti; ) exp [ — i 52 o-i 



The effective Planck constant is hcs — T — 2tt/N] Si = ujiT, S2 — W2T; e is the coupling 
strength between the qubits and the environment. The classical limit HeS ^ is obtained 
by taking T ^ and k -^ 00, in such a way that K = kT is kept fixed. 

I am interested in the case in which the environment (the kicked rotator) is chaotic 
(that is, with K ^ 1). The two qubits arc initially prepared in a maximally entangled 
state, so that they are disentangled from the environment. Namely, I suppose that at 
t = the system is in the state 

(72) 1*0) = 1'/'+) ® Ho) , 

where \(j)'^) = -^ (|00) + |11)) is a Bell state (the particular choice of the initial maximally 
entangled state is not crucial for what follows), and Jt/iq) = X^n Cn\n) is a generic state of 
the kicked rotator, with c„ random coefficients such that '^^ |c„p = 1, and \n) eigenstates 
of the momentum operator. The evolution in time of the global system (kicked rotator 
plus qubits) is described by the unitary operator U defined in Eq. (71). Therefore, any 
initial pure state |^o) evolves into another pure state |^(i)) = t/*|*o)- The reduced 
density matrix pi2{t) describing the two qubits at time t is then obtained after tracing 
|\l'(t))(^(i)| over the kicked rotator's degrees of freedom. 

In the following I will focus my attention on the time evolution of the entanglement of 
formation (see Appendix D) E12 between the two qubits and that between them and the 
kicked rotator, measured by the reduced von Neumann entropy 6*12 = — Tr [pi2 logj P12] 
of the reduced density matrix pi2. Clearly, for states like the one in Eq. (72), we have 
-£'12(0) = 1, S'i2(0) = 0. As the total system evolves, we expect that E12 decreases, while 
'S'12 grows up, thus meaning that the two-qubit system is progressively losing coherence. 

If the kicked rotator is in the chaotic regime and in the semiclassical region hcn <C 1 , 
it is possible to drastically simplify the description of the system in Eq. (66) by using 
the random phase-kick approximation, in the framework of the Kraus representation 
formalism. Since, to a first approximation, the phases between two consecutive kicks 
in the chaotic regime can be considered as uncorrelated, the interaction with the envi- 
ronment can be simply modeled as a phase-kick rotating both qubits through the same 
random angle about the z-axis of the Bloch sphere. This rotation is described in the 
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Fig. 6. - Reduced von Neumann entropy 512 (main figure) and entanglement E12 (inset) as a 
function of time at K ^ 99.73, 5i — 10^^, 82 = a/2(5i, e — 8x10"^. The thin curves correspond to 
different number of levels for the environment (the kicked rotator) (iV = 2^, 2^", 2", 2^^, 2", 2" 
from bottom to top in the main figure and vice versa in the inset). The thick curves give the 
numerical results from the random phase model (74). 



{|00), |01), |10), |11)} basis by the unitary matrix 
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where the angle 6 is drawn from a uniform random distribution in [0, 2tt). The one- kick 
evolution of the reduced density matrix pi2 is then obtained after averaging over 9: 



(74) 



P12 



1 

2^ 



deR(e) e-'^'"^^"^ e-'^'"^^'^ P12 e'^"'^^"e'^^''^^'^R'i{e). 



In order to assess the validity of the random phase-kick approximation, model (66) 
is numerically investigated in the classically chaotic regime K :^ 1 and in the region 
heS <C 1 in which the environment is a semiclassical object. Under these conditions, we 
expect that the time evolution of the entanglement can be accurately predicted by the 
random phase model. Such expectation is confirmed by the numerical data shown in 
Fig. 6. Even though differences between the two models remain at long times due to 
the finite number N of levels in the kicked rotator, such differences appear at later and 
later times when N ^ 00 (ficff ^ 0). The parameter K has been chosen much greater 
than one, so that the classical phase space of the kicked rotator can be considered as 
completely chaotic. Note that the value K « 99.72676 is chosen to completely wipe off 
memory effects between consecutive and next-consecutive kicks (see Ref. [72] for details) . 
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I point out that the random phase model can be derived from the Caldeira-Leggett 
model with a pure dephasing coupling ex {ai +ai ) ^j, gkQk, with gi~ coupling constant 
to the fc-th oscillator of the environment, whose coordinate operator is qk [76, 81]. This 
establishes a direct link between the chaotic single-particle environment considered in 
this paper and a standard many-body environment. 

10. Final remarks 

The role of entanglement as a resource in quantum information has stimulated in- 
tensive research aimed at unveiling both its qualitative and quantitative aspects. The 
interest is first of all motivated by experimental implementations of quantum informa- 
tion protocols. Decoherence, which can be considered as the ultimate obstacle in the 
way of actual implementation of any quantum computation or communication protocol, 
is due to the entanglement between the quantum hardware and the environment. The 
decoherence-control issue is expected to be particularly relevant when the state of the 
quantum system is complex, namely when it is characterized by a large amount of multi- 
partite entanglement. It is therefore important, for applications but also in its own right, 
to scrutinize the robustness and the multipartite features of relevant classes of entangled 
states. In this context, random states play an important role, both for applications in 
quantum protocols and in view of a, highly desirable, statistical theory of entanglement. 

Such studies have deep links with the physics of complex systems. In classical physics, 
a well defined notion of complexity, based on the exponential instability of chaos, exists, 
and has profound links with the notion of algorithmic complexity [82, 83]: in terms of the 
symbolic dynamical description, almost all orbits are random and unpredictable. On the 
other hand, in spite of many efforts (see [84] and references therein) the transfer of these 
concepts to quantum mechanics still remains elusive. However, there is strong numerical 
evidence that quantum motion is characterized by a greater degree of stability than 
classical motion (see [85, 36, 7, 86]). This has important consequences on the stability 
of quantum algorithms; for instance, the robustness of the multipartite entanglement 
generated by chaotic maps and discussed in Sec. 7 is related to the power-law decay 
of the fidelity time scales for quantum algorithms [87, 88, 54, 89] which, in turn, is a 
consequence of the discreteness of the phase space in quantum mechanics [85, 36, 7, 86]. 
If we consider the chaotic classical motion (governed by the Liouville equation) of some 
phase-space density, smaller and smaller scales are explored exponentially fast. These 
fine details of the density distribution are rapidly lost under small perturbations. In 
quantum mechanics, there is a lower limit to this process, set by the size of the Planck 
cell, and this reduces the complexity of quantum motion as compared to classical motion. 

Finally, the fundamental, purely quantum notion of entanglement is expected to play 
a crucial role in characterizing the complexity of a quantum system [90] . I believe that 
studies of complexity and multipartite entanglement will shed some light on series of 
very important issues in quantum computation and in critical phenomena of quantum 
many-body condensed matter physics. 

Appendix A. 

Separability criteria 

A'l. The Peres criterion. - The Peres criterion [91] provides a necessary condition 
for the existence of decomposition (17), in other words, a violation of this criterion is 
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a sufficient condition for entanglement. This criterion is based on the partial transpose 
operation. Introducing an orthonornial basis {|i)A|Q^)B} in the Hilbert space Ti.AB £^s- 
sociated with the bipartite system A + B, the density matrix pab has matrix elements 
{pAB)ia;ji3 = A (*|s (ct|pAB | j) A |/?)b • The partial transpose density matrix is constructed 
by only taking the transpose in either the Latin or Greek indices (here Latin indices refer 
to Alice's subsystem and Greek indices to Bob's). For instance, the partial transpose 
with respect to Alice is given by 

Since a separable state pab can always be written in the form (17) and the density 
matrices pAk and p^fe have non-negative eigenvalues, then the overall density matrix 
PAB also has non- negative eigenvalues. The partial transpose of a separable state reads 

(A.2) P^A%'^^PkP^Ak®PBk- 

k 

Since the transpose matrices p^^, = p^j, are Hermitian non-negative matrices with unit 
trace, they are also legitimate density matrices for Alice. It follows that none of the 
eigenvalues of p^-^ is non-negative. This is a necessary condition for decomposition (17) 
to hold. It is then sufhcient to have at least one negative eigenvalue of p^-^ to conclude 
that the state pab is entangled. 

It can be shown [92] that for composite states of dimension 2x2 and 2x3, the Peres 
criterion provides a necessary and sufficient condition for separability; that is, the state 
Pab is separable if and only if p^g is non-negative. However, for higher dimensional 
systems, states exist for which all eigenvalues of the partial transpose are non-negative, 
but that are non-separable [93]. These states are known as bound entangled states since 
they cannot be distilled by means of local operations and classical communication to 
form a maximally entangled state [94] . 

I stress that the Peres criterion is more sensitive than Bell's inequality for detecting 
quantum entanglement; that is, there are states detected as entangled by the Peres 
criterion that do not violate Bell's inequalities ([91]). 

A'2. Entanglement witnesses. - A convenient way to detect entanglement is to use 
the so-called entanglement witnesses [92, 95]. By definition, an entanglement witness is 

a Hermitian operator W such that 1y:(W p^g^ > for all separable states Pab ''^t'.ilc 

there exists at least one state Pab such that Tr(T4^p^^^) < 0. Therefore, the negative 

expectation value of T4^ is a signature of entanglement and the state p^g is said to be 
detected as entangled by the witness W. 

The existence of entanglement witnesses is a consequence of the Hahn-Banach theo- 
rem: Given a convex and compact set S and pab ^ S, there exists a hyperplane that 
separates pab from 5*. 

This fact is illustrated in Fig. 7. The set S of separable states is a subset of the set of 
all possible density matrices for a given system. The dashed line represents a hyperplane 
separating an entangle state pab from S. The optimized witness Wopt (represented by 
a full line) is obtained after performing a parallel transport of the above hyperplane, so 
that it becomes tangent to the set of separable states. Therefore, the optimized witness 
Wopt detects more entangled states than before parallel transport. Note that, in order to 
fully characterize the set S of separable states one should find all the witnesses tangent 
toS". 



30 



G. BENENTI 




Fig. 7. - Schematic drawing of entanglement witnesses. 



The concept of entanglement witness is close to experimental implementations and 
detection of entanglement by means of entanglement witnesses has been realized in several 
experiments [96, 97, 98]. The more negative expectation value of entanglement witness 
we find, the easier it is to detect entanglement of such a state. The expectation value 
of W also provides lower bounds to various entanglement measures [99, 100]. Finally, 
it is interesting to note that violation of Bell's inequalities can be rewritten in terms of 
non-optimal entanglement witnesses [95, 101]. 

In general, classification of entanglement witnesses is a hard problem. However, much 
simpler is the issue with the so-called decomposable entanglement witnesses [102]. By 
definition, a witness is called decomposable if 



(A.3) 



W^P + Q^ 



P, Q > 0, 



that is, with positive semidefinite operators P, Q. Decomposable entanglement witnesses 
can only detect entangled states with at least one negative eigenvalues of p^^ . Therefore, 
similarly to the Peres criterion, decomposable witnesses do not detect bound entangled 
states. Note, however, that entanglement witnesses are closer to experimental implemen- 
tations than the Peres criterion, as full tomographic knowledge about the state is not 
needed. 
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Purity of random states 

Let us write a iV-level random state in the form 



N-l 



(B.l) 



|V^) = ^rfee'^'=|fc), 



k=0 
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where (j)k are independent random variables uniformly distributed in [0, 27r) and r = 
(rg, ...,rjv-i) is a random point uniformly distributed on the unit hypersphere §^^^ = 
{r e M^|r^ = 1}, with distribution function 

AT-l 

(B.2) p(r) = Cat n ru5{v^ - 1), 

fc=0 

with Cn normalization constant to be determined later. 

Given a bipartition of the Hilbert space of the system into two parts, A and B, with 
dimensions Na and Nb, the purity reads 

Na-INb-I 

(B.3) P= ^ Yl '^Ji'^J'i'^J'i'^31' exp[i((/'j7 - 0// + (t)j'v - (t)ji')], 

j.j'=0 l,l'=0 

where |fc) = \jl) — \j)a ® |0s- Following [43], I split P in two parts: 

(B.4) P^X + M, 

where 

(B.5) ^ = X! ' X! ''^Ji^3'i'rj'i''rji' exp[«(0j7 - (t)j'i + (t^j'v - (t)jv)\, 

3,3' 1,1' 



(B.6) m^j:'T. 44i + T.T. '44' + E 

3,3' I 3 1,1' 3,1 



^3l 



where ^ means that equal indexes are banned in the sum. Since {e'"^'') = we obtain 
{X) =. 0. Therefore, 

(B.7) P = (M) = N{Na + ^b - 2)(r2r2> + N{tI), 

where I have used (r^) — {r^) for all k and {r^r'j.,) — {r^rl) for all k, k' with k ^ k' . 
I now evaluate the marginal distribution 

p{ro,...,r,n-i) ^ CnTo- --rm-i drmTm drm+iTm+i--- drN-irN-i5{r'^ -I) 

Jo Jo Jo 



-TT i~o--- r,n-i / drmr„i ■■ drM-2rN-2 

^ Jo Jo 



. Af-m-1 



32 G. BENENTI 

In particular, 

The normalization condition L drop{ro) = 1 allows us to determine Cat = 2^ {N — 1)!. 
Thus, we obtain 



(B.IO) p{ro) ^ 2{N ^ l)ro{l - r'o) 



2\N-2 



(B.ll) (rl) = / drorlpiro) = 



" ' NiN + iy 



(B.12) 



p{ro,ri) - 4{N - 1){N - 2)rori(l - r^ - rjf-^, 



(B.13) ^''^'■'^"i '^'^°''°i ^^i'^iP(^o,ri)^ ^^^^^^^ . 

After substitution of Eqs. (B.ll) and (B.13) into (B.7) we readily obtain Lubkin's formula 
(23). 

The variance ap can be computed with the same technique as above [43, 35]. However, 
to obtain the variance (24) for large N it is suifcient to replace in Eqs. (B.5) and (B.6) 
rfe with its mean value l/y/N: 

(B.14) 4 = (P2) -Pl = {X^) + (M^) - P| « {X^) « A. 

We can see from Eqs. (B.5) and (B.6) that X and M are sums of 0{N'^) terms of order 
1/N'^. Therefore, the central limit theorem implies that, for large iV, the purity tends to 
a Gaussian distribution with mean Pl and variance ap. Finally, I note that all moments 
of the purity have been recently computed [103]. 
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The sawtooth map 

The sawtooth map is a prototype model in the studies of classical and quantum 
dynamical systems and exhibits a rich variety of interesting physical phenomena, from 
complete chaos to complete integrability, normal and anomalous diffusion, dynamical 
localization, and cantori localization. Furthermore, the sawtooth map gives a good ap- 
proximation to the motion of a particle bouncing inside a stadium billiard (which is a 
well-known model of classical and quantum chaos) . 
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C'l. Classical dynamics. - The sawtooth map belongs to the class of periodically 
driven dynamical systems, governed by the Hamiltonian 

2 +°o 

(c.i) H{e, n-T) = '^ + v{e) J2 s(^ - JT) , 

j = -oo 

where {n,9) are conjugate action-angle variables {0 < < 27r). This Hamiltonian is the 
sum of two terms, H{9, n: r) = Ho{n) + U{6; t), where Ho{n) = n^/2 is just the kinetic 
energy of a free rotator (a particle moving on a circle parametrized by the coordinate 6), 
while 

(C.2) U{9-r)^V{0)Y,Hr-jT) 

j 

represents a force acting on the particle that is switched on and off instantaneously at 
time intervals T. Therefore, we say that the dynamics described by Hamiltonian (C.I) 
is kicked. The corresponding Hamiltonian equations of motion are 

an dV{0) S 

. dH 

1 = — — = n. 
on 

These equations can be easily integrated and one finds that the evolution from time tT~ 
(prior to the t-th kick) to time (t+ 1)T^ (prior to the (t+ l)-th kick) is described by the 
map 

r fit+i ^nt+ F{0) , 
(C.4) 

(^ 6*4+1 =dt+ Tfit+i , 

where the discrete time t — t jT measures the number of map iterations and F(ff) — 
—dV{6)/d6 is the force acting on the particle. 

In the following, we focus on the special case V{6) ~ —k{9 — 7r)^/2. This map is 
called the sawtooth map, since the force F(9) = —dV{9)/d9 = k{9 — n) has a sawtooth 
shape, with a discontinuity at ^ = 0. For such a discontinuous map the conditions of 
the Kolmogorov-Arnold-Moser (KAM) theorem are not satisfied and, for any fc ^ 0, the 
motion is not bounded by KAM tori. By rescaling n -^ I = Tn, the classical dynamics 
is seen to depend only on the parameter K — kT. Indeed, in terms of the variables (/, 9) 
map (C.4) becomes 

r It+i ^It+ K{9 - n) , 
(C.5) 

I 0t+i ~ 9t + It+i ■ 

The sawtooth map exhibits sensitive dependence on initial conditions, which is the dis- 
tinctive feature of classical chaos: any small error is amplified exponentially in time. In 
other words, two nearby trajectories separate exponentially, with a rate given by the 
maximum Lyapunov exponent A, defined as 

(C.6) A= lim -Inf^ 

|t|^oo t \0o 
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where 5t = ^/ISIt]-^ + [Sdt]'^. To compute Sit and 66t, we differentiate map (C.5), obtain- 
ing 



(C.7) 



SIt+i 
59t+i 



M 



Sit 

set 



1 K 

1 1 + K 



Sit 

set 



The iteration of map (C.7) gives Sit and Set as a function of SIq and Seo {SIq and SOq 
represent a change of the initial conditions). The stabihty matrix M has eigenvalues 
^{2 + K zL y/lO+4K), which do not depend on the coordinates / and e and are 
-4 < X < and real for K < -A and K > 0. Thus, the classical 
4 < i^ < and completely chaotic for K < —4 and K > 0. 
.)* asymptotycally in t, and therefore the maximum Lyapunov 
Similarly, we obtain A = In |/i_| for K < —4. In the stable region 



complex conjugate for 
motion is stable for - 
For K > 0, St oc (^_ 
exponent is A = In /i_|_ . 



-4<X<0, A = 0. 

The sawtooth map can be studied on the cylinder [/ G (— oo, +oo)], or on a torus of 
sinite size (— ttL < I < ttL, where L is an integer, to assure that no discontinuities are 
introduced in the second equation of (C.5) when / is taken modulus 2'kL). Although the 
sawtooth map is a deterministic system, for iiT > and K < —4 the motion along the 
momentum direction is in practice indistinguishable from a random walk. Thus, one has 
normal diffusion in the action (momentum) variable and the evolution of the distribution 
function /(/, t) is governed by a Fokker-Planck equation: 



(C.8) 



dl^d_ /I a/ 

dt dl\2 dl 



The diffusion coefficient D is defined by 



(C.9) 



D = lini 

t — )-00 



t 



where A/ = I — {I), and (. . .) denotes the average over an ensemble of trajectories. If 
at time t — we take a phase space distribution with initial momentum Iq and random 
phases < 6 < 2tt, then the solution of the Fokker-Planck equation (C.8) is given by 



(CIO) 



f{I,t)^ 



1 



V27rDi 



cxp 



(I-Io? 



2Dt 



The width ■^/((A/f)^) of this Gaussian distribution grows in time, according to 



(C.ll) 



((A/*: 



D{K)t. 



For K > 1^ the diffusion coefficient is well approximated by the random phase approxi- 
mation, in which we assume that there are no correlations between the angles (phases) 
6 at different times. Hence, we have 



{G.12)D{K) w ((A/(i)2 



1 

2^ 



de{M 



2ti 7o 3 



K' 



where A/^^^ — It+i — It is the change in action after a single map step. For Q < K < 1 
diffusion is slowed, due to the sticking of trajectories close to broken tori (known as 
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cantori), and we have D{K) « 3.3 K^^^ (this regime is discussed in Ref. [104]). For 
—4 < K < the motion is stable, the phase space has a complex structure of elliptic 
islands down to smaller and smaller scales, and one can observe anomalous diffusion, that 
is, ((AJ)^) ex i", with a 7^ 1 (see Ref. [54]). The cases K — —1, —2, --3 are integrable. 

C'2. Quantum dynamics. - The quantum version of the sawtooth map is obtained by 
means of the usual quantization rules, — > and n ^ h = —id/dO (we set h — 1). The 
quantum evolution in one map iteration is described by a unitary operator U , called the 
Floquet operator, acting on the wave vector \ip): 



(C.13) lV>)t+i = [> IV-)* = exp 



M+i)T- 
i / dTH{0,i;T) 

JIT- 



IV-)*, 



where H is Hamiltonian (C.l). Since the potential V{d) is switched on only at discrete 
times IT, it is straightforward to obtain 

(C.14) lV')m-e-^"'/'e-''^''^l^)„ 

which for the sawtooth map is just Eq. (46). It is important to emphasize that, while 
the classical sawtooth map depends only on the rescaled parameter K = kT, the corre- 
sponding quantum evolution (C.14) depends on k and T separately. The effective Planck 
constant is given by fi-cff = T. Indeed, if we consider the operator / = Tn (/ is the 
quantization of the classical rescaled action /), we have 



(C.15) [9, i] = T[§, n]::^iT ^ ih, 



-off- 



The classical limit ftcff ^ is obtained by taking k ^ 00 and T ^ 0, while keeping 
K = kT constant. 

In the quantum sawtooth map model one can observe important physical phenomena 
like dynamical localization [105]. Indeed, due to quantum interference effects, the chaotic 
diffusion in momentum is suppressed, leading to exponentially localized wave functions. 
This phenomenon was first found and studied in the quantum kicked-rotator model [106] 
and has profound analogies with Anderson localization of electronic transport in disor- 
dered materials [107]. Dynamical localization has been observed experimentally in the 
microwave ionization of Rydberg atoms [108] and in experiments with cold atoms [109]. 
In the quantum sawtooth map also cantori localization takes place: In the vicinity of a 
broken KAM torus, a cantorus starts to act as a perfect barrier to quantum wave packet 
evolution, if the flux through it becomes less than h [110, 111, 112, 113, 114]. 

C'3. Quantum algorithm. - In the following, we describe an exponentially efficient 
quantum algorithm for simulation of the map (46) [54]. It is based on the forward/ 
backward quantum Fourier transform between momentum and angle bases. Such an 
approach is convenient since the Floquet operator [/, introduced in Eq. (46), is the 
product of two operators, IJu = e^i^is-^) /2 ^j^^j jj^ — g-«^" /2^ diagonal in the 9 and 
n representations, respectively. This quantum algorithm requires the following steps for 
one map iteration: 

• Apply Uk to the wave function ip{6). In order to decompose the operator Uk into 



36 G. BENENTI 

one- and two-qubit gates, we first of all write 9 in binary notation: 
(C.16) e = 2TT^aj2-\ 

with Ui G {0,1}. Here Uq is the number of qubits, so that the total number of 
levels in the quantum sawtooth map is A'^ = 2""' . From this expansion, we obtain 

^ ' ^ ' ^ \2i^ 2nJ \2i^ 2n 

This term can be put into the unitary operator Uk, giving the decomposition 



(C.18) ^^kie-.f/2^ -Q c^p 



2J1 2nJ\ 2J2 2n 



which is the product of n^ two-qubit gates (controlled phase-shift gates), each acting 
non-trivially only on the qubits ji and j2. In the computational basis {{ctj-^aj^} — 
|00), |01), |10), |11)} each two-qubit gate can be written as cxp{i2'ii'^kDj^,j^)^ where 
^3i,J2 is a diagonal matrix: 
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• The change from the 9 to the n representation is obtained by means of the quantum 
Fourier transform, which requires Uq (single-qubit) Hadamard gates and ■2nq{nq^l) 
(two-qubit) controlled phase-shift gates (see, e.g., [7, 8]). 

• In the n representation, the operator Ut has essentially the same form as the oper- 
ator Uk in the 9 representation and can therefore be decomposed into n^ controlled 
phase-shift gates, similarly to Eq. (C.18). 

• Return to the initial 9 representation by application of the inverse quantum Fourier 
transform. 

Thus, overall, this quantum algorithm requires '^rii + Uq gates per map iteration (3n^ — n^ 
controlled phase-shifts and 2nq Hadamard gates). This number is to be compared with 
the 0(nq2"'') operations required by a classical computer to simulate one map iteration 
by means of a fast Fourier transform. Thus, the quantum simulation of the quantum 
sawtooth map dynamics is exponentially faster than any known classical algorithm. Note 
that the resources required to the quantum computer to simulate the evolution of the 
sawtooth map are only logarithmic in the system size N . Of course, there remains the 
problem of extracting useful information from the quantum computer wave function. For 
a discussion of this problem, see Refs. [7, 115]. Finally, I point out that the quantum saw- 
tooth map has been recently implemented on a three-qubit nuclear magnetic resonance 
(NMR)-based quantum processor [116]. 
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Appendix D. 

Entanglement of formation and concurrence 

Any state p can be decomposed as a convex combination of projectors onto pure 
states: 

(D.i) p = Ep^i^^)(^^-I- 

3 

Tire entanglement of formation E is defined as the mean entanglement of the pure states 
forming p, minimized over all possible decompositions: 

(D.2) E{ps)^MTp,E{\^,)), 



dec 



where the (bipartite) entanglement of the pure states \tpj) is measured according to 
Eq. (9). 

The entanglement of formation £^12 of a generic two-qubit state pi2 can be evaluated 
in a closed form following Ref. [117]. First of all we compute the concurrence, defined 
as C = max(Ai — A2 — A3 — A4, 0), where the Aj's are the square roots of the eigenvalues 
of the matrix R — pi2pi2, in decreasing order. Here pi2 is the spin flipped matrix of 
P12, and it is defined by pi2 = {<Jy ® cry) Pi2 (cFy (E" cFy) (note that the complex conjugate 
is taken in the computational basis {|00), |01), |10), |11)}). Once the concurrence has 
been computed, the entanglement of formation is obtained as _E == h{{l + ^/l~~C^)/2), 
where h is the binary entropy function: h{x) — — xloggX — (1 — a:;)log2(l — x\ with 
a; = (1 + Vl-C2)/2. 

The concurrence is widely investigated in condensed matter physics, in relation to the 
general problem of the behavior of entanglement across quantum phase transitions [118]. 
For studies of the relation between entanglement and integrability to chaos crossover in 
quantum spin chain, see [119, 120, 121, 122, 123, 84], and references therein. 
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